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Abstract. Tunnelling processes through black hole horizons have recently been investigated in the frame- 
work of WKB theory discovering interesting interplay with the Hawking radiation. A more precise and 
general account of that phenomenon has been subsequently given within the framework of QFT in curved 
spacetime by two of the authors of the present paper. In particular, it has been shown that, in the limit 
of sharp localization on opposite sides of a Killing horizon, the quantum correlation functions of a scalar 
field appear to have thermal nature, and the tunnelling probability is proportional to exp{— finawkingE}. 
This local result is valid in every spacetime including a local Killing horizon, no field equation is nec- 
essary, while a suitable choice for the quantum state is relevant. Indeed, the two point function has to 
verify a short-distance condition weaker than the Hadamard one. In this paper we consider a massive 
scalar quantum field with a (f> 3 self-interaction and we investigate the issue whether or not the black hole 
radiation can be handled at perturbative level, including the renormalisation contributions. We prove 
that, for the simplest model of the Killing horizon generated by the boost in Minkowski spacetime, and 
referring to Minkowski vacuum, the tunnelling probability in the limit of sharp localization on opposite 
sides of the horizon preserves the thermal form proportional to exp{— /3hE} even taking the one- loop 
renormalisation corrections into account. A similar result is expected to hold for the Unruh state in the 
Kruskal manifold, since that state is Hadamard and looks like Minkowski vacuum close to the horizon. 



1 Introduction 

1.1. Hawking radiation as tunnelling process. In remarkable papers Parikh and Wilczek 
|PW00] and, independently, Volovik [Vo99] computed the transition probability Pe, also called 
tunnelling probability, of a particle with energy E through a Schwarzschild Black Hole horizon. 
They found the very thermal behavior Q with the characteristic Hawking temperature Tjj (it 
is throughout assumed kBoitzmann = 1)- This brilliant result implies several physically notable 
consequences and proposes a new viewpoint on the black hole radiation |Ha75| issue. In par- 
ticular, since no detection of radiation at future null infinity is exploited, the outcome of the 
computation suggests that the black hole radiation can be viewed as a completely local phe- 
nomenon. Furthermore, just what happens in a neighborhood of a point on the horizon seems 



to actually matter and thus no global black hole structure seems to be necessary in this picture. 
The mentioned tunnelling probability: 

?E = \ {lpEx\lpEy)\ 2 

is computed for one particle with energy E, where the wavefunctions ipEx, ^Ey are localized 
around x = (t,ri,0,<p) and y = (t,r2,9,(p) respectively and these events are separated by the 
horizon. The understood quantization procedure is performed referring to the Painleve time 
t appearing in the corresponding explicit expression of the Schwarzschild metric. The overall 
authors' idea is to take advantage of WKB method to approximate Ye and to study the leading 
order approximation for the case of x infinitesimally close to y, though separated from it by the 
horizon: 

= lim \ {ipEx\ipEy)\ 2 ~ hm 

y— >x y^rx 

However a snag pops up: The integral in the exponent diverges as soon as x approaches y. In 
|PW00] this difficulty is turned out into an advantage by exploiting a suitable complex plane 
Feynman-like regularization. In this way an imaginary part arises in the integral yielding: 

This, nowadays very popular, result has been subsequently reproduced by various authors clean- 
ing some unclear technical issues [APS08, APGS08J but also generalizing the geometrical outset, 
even quite considerably, encompassing new physically remarkable situations like dynamical black 
holes horizons and referring to other kinds of particles and in relation to other physical issues 
as the back reaction on the metric |CV991 IANVZ051 IMV051 IDNVZZ071 IKM071 IHDVNZ09] (see 
[VADTT] for a survey). However, the presented machinery seems to remain plagued by some 
unresolved problems analysed in [MP 12]. First of all, the appearance of Tjj seems to be suspi- 
ciously related with the choice of the adopted complex-plane regularization procedure. Though, 
on the other hand, the correct value of the temperature very hardly shows up by chance. Fur- 
thermore some key tools as the notions of particle (but also time and energy) are ambiguously 
defined in curved spacetime, due to the absence of the Poincare symmetry. Finally, despite it 
is clearly suggested by the flavour of the final result, it is by no means clear how the result is 
independent from the whole Black Hole geometry. This is because was obtained in |PW00] 
dealing with the Schwarzschild black hole manifold. 

1.2. The viewpoint of algebraic QFT in curved spacetime. The proper framework of 
algebraic QFT in curved spacetime was adopted in |MP12| to clarify the physical meaning of 
Parikh-Wilczek's result. Let us review the outcome of that analysis, addressing to [Wa84] for all 
geometric notions we employ. In a 4-dimensional smooth spacetime M with Lorentzian metric 
g having signature —,+,+,+, we consider an open set, C M, where a smooth vector field K 
exists satisfying the following requirements. 
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(a) K is a Killing field for g in 0. 

(b) contains the local Killing horizon IK i.e. a 3-submanifold invariant under the action 
of K with K a K a = on IK. 

(c) The orbits of X in are diffeomorphic to an open interval / and topologically IK = Jx 2 
(23 being a 2-dimensional cross section). 

(d) The surface gravity k ^ is constant on IK. (k is defined by V a (Ki ) K b ) = —2nK a .) 

We shall make use of a standard null coordinate system U, V, s adapted to IK, where U G I is the 
affine parameter of the null geodesies forming K, y is the affine parameter of the null geodesies 
crossing K once - with the choice of the origin such that x G K iff V(x) = - and s denotes a 
pair of coordinates over IB where U = 0. We refer to [MP12J for a precise definition. 
As the computation will not depend on the geometry outside 0, the horizon may (smoothly) 
cease to exist outside 0. The requirement k = constant along K means that the thermodynamic 
equilibrium has been reached on K, since k = 2ttTh- Notice that conditions (a)-(d) are quite 
general. For example they are satisfied around points of the future horizon of a non-extremal 
black hole in the Kerr-Newman family, including charged and rotating black holes. Furthermore, 
they are also valid both for "realistic" black holes produced by collapsed matter - so that only 
the future horizon exists - and even for eternal black holes - whose manifolds include white hole 
regions as in Kruskal spacetime. Finally, our picture includes also situations where the collapse 
starts, reaches a sort of local equilibrium and it stops after a while, without giving rise to a 
complete BH structure. 

Having discusses the geometric setup we pass now to present quantum matter we consider. 
From now on A is the unital *-algebra generated by abstract scalar field operators 4>(f) 
with / G Co° (M) (the space of smooth complex and compactly supported functions on M) such 
that: 

(Rl) cf>(af + bf')_= acf>(f) + bct>(f) if a, b G C and /, /' G C °°(M); 
(R2) «/>(/)* = </>(/) for / G C^(M); 

(R3) [<fi(f), 4>(f')] = f° r causally disjoint supp(f), supp(f'). 

Notice that, among these requirements no field equation is assumed. However, since we intend 
to compute the correlation function o;(^(/)(/>(/ / )) with respect to the state uj on A, when 
supp(f),supp(f) C are "very close" to K, a condition on the state is necessary. 
Supposing that is geodesically convex, the signed squared geodesic distance of x and y, 
o~(x,y), is thereon well-defioned, t is any timelike coordinate and a standard e — > + prescription 
is assumed whenever indicated. With these notations, we assume that the integral kernel u>2 of 
uj((fi(f)(j)(f')) verifies the last hypothesis we need. 

(R4) The short-distance behaviour holds: 

/ /\ D(p,p') 

U2{P,P) ■= —, 77 +W e (p,p) 
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where a e (p,p') := o-(p,p') + 2ie(t p — t p /) + e 2 , D is smooth and admits a constant c / 
with D(p,p') = c if p, p' G !K and o~(p,p') = Finally 7u e is a distribution less singular 
than l/cr e . 

As in [MP 12] we say that w t is less singular than l/<7 e if for every e, w; e are measurable functions 
such that: 

(i) w t (p,p') — > w'(p,p'), almost everywhere in (p,p r ) as e — >• + for some function w' and w e is 
e-uniformly bounded by a locally M 2 -integrable function; 

(ii) v/(V, U, s, V , U', s') -> u/'(£7, s, £/', s') almost everywhere in ([/, s, s') when (V, V) -> 
(0, 0) for some function w" on !K 2 that is (V, V' )-uniformly bounded by a locally CK 2 -integrable 
function. 

Remark 1.1. (1) An example of w e less singular than l/cr e is, for every fixed 5 > 0: 



w e = t~¥~4^& + fe(x,y) In a e (x,y) + r e (x,y) (2) 



where for fixed e 6 (0, eo), f e , r e , h e are measurable with \f e (x,y)\ < K \g t (x,y)\ < K' 
and \h e (x, y)\ < K" for constants K,K',K' and all (e,x,y) £ (0, eo) x x and, finally, 
f e (x,y) -)■ f(x,y), r e (x,y) -)■ r(x,y) and h e (x,y) -)• for e ->• 0+ if cr(x,y) / 0. Above 

and throughout the cuts in the complex domain of In z and z Q with a Z are supposed to stay 
along the negative real axis. 

We stress that Gaussian Hadamard states for a linear Klein-Gordon field eft satisfy the require- 
ment in (R4) with w e as in ^ with h € = 0. 

(2) A requirement like (R4) was essentially assumed in |FH90| . dealing with a free scalar field, 
to prove that this class of states produces the black hole radiation at future null infinity for a 
spacetime containing spherically symmetric collapsing matter giving rise to a black hole. 

(3) If we assume that oj is Gaussian and supp(f), supp(f') are separated by the horizon, then 
in the Hilbert space of the GNS representation of u, \uj{(j)(f)(f>{f'))\ 2 is nothing but the transi- 
tion probability of a pair of one-particle states localized on the opposite sides of the horizon: a 
tunnelling probability. 

We re- interpret the limit "x — > y" in Parikh-Wilczek picture as follows: 

lim ui^UMfx)) = Km u; 2 (f x , f x ) (3) 

A— >Q+ A— s>0+ 

where, as usual, on the right-hand side we denoted the distribution with the same symbol of its 
integral kernel. Furthermore, the limit A — > + shrinks the supports of f and f on "K. Explicitly, 
making use of the previously mentioned null coordinate system U, V, s adapted to "K: 

f x (V,U,x) = \f(j,U,s) . (4) 



1 Apparently, a stronger requirement on D shows up in |MP12j . actually this only requirement was exploited. 
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To remove an infrared divergence arising in the computation of co((j)(f\)(j)(f' :x )) as A — >■ + , we 
also assume that: 

/ = d v F and /' = d v F' for a couple F, F' G C§°(0) (5) 

Finally we need to specify the notions of time and energy, for the (locally) external region at 
least where K is timelike. Therein V ~ e~ KT . The natural choice for time is the the parameter 
t of the Killing field K. The (locally) internal region is not stationary (K is spacelike and 
V ~ —e KT ) so no geometrically natural notion of time can exists there. Therefore we stick with 
r in that region, too. The Fourier transform with respect to r defines the energy spectrum with 
respect to the notion of energy E associated with r. We can state the main result of |MP12j 
that also includes the case of test functions both supported in the external region. 

Theorem 1.1. Assuming that for the open set OcM the hypotheses (a)-(d) hold and, referring 
to a state oj on A, that the requirements (R1)-(R4) hold true as well, one has: 



A hm W (0(AM/ A )) = — J J ginh(/W2) EdEdUdUd^x) 



where f and f'have supports separated by the horizon, verify j^J), and F(E, U, x), F'(E, U, x) 
denotes the r-Fourier transform of F(e~ KT ,U, x) and F'{— e~ KT , U, x) respectively. For wave 
packets sharply centered around a large value Eq of the energy, ( flip yields: 

lim |u^(/a)<M/a))| 2 ~ CEl e~^ Eo . 

A— >0 

Finally, if both the supports of f and f stay in the region where the Killing vector is timelike, 
a proper Bose spectrum arises: 



hm + u]{(/)(f x )(t>(f x )) = — / / _ ^_ PhE EdEdUdU dn{x) 



In this way, adopting the viewpoint of algebraic QFT in curved spacetime, Parikh and Wilczek's 
result acquires a precise and rigorous meaning. As our computation is completely local in space 
and time, it strongly supports the idea that the Hawking radiation is ( also ) a local phenomenon, 
independent from the existence of a whole black hole^j The result is independent form the state 
of the quantum field, provided it belongs to a large class including the physically significant 
Hadamard states. That class of states enjoys physically fundamental properties in developing 
QFT in curved spacetime, actually in a semiclassical treatment of quantum gravity |Wa94) 
the back-reaction on the metric can be computed because they admit a well-defined stress 
energy tensor jMo03] . moreover, the renormalisation procedure can be successfully implemented 
referring to them jBFOOl IHoWaOH IHoWa021 lBDF09j . 



2 However, as stressed in [MP 12] without fixing the value of K at some point, for instance at infinity, a constant 
rescaling ambiguity remains in the definition of K affecting the value of Th- In an asymptotically flat spacetime 
this ambiguity is removed assuming that K defines the Minkowski standard time far away from the horizon. 
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2 The simplest interacting case: £/ = in Rindler spacetime 

2.1. Idea and motivation. As the result of |MP12] does not depend on any field equation, 
but only on the short distance behavior (R4) of the two-point function of the considered state, it 
must be true even considering interacting quantum fields provided a suitable scaling limit of the 
two-point function holds [FH87, Bu96]. Even in that case, it is however far from obvious that 
(R4) also holds when treating QFT perturbatively, taking the renormalisation corrections into 
account. As a matter of fact we will focus on the simplest model given by the Lj = Jy(^> 3 self- 
interaction in Minkowski spacetime, referring to the Killing horizon generated by a boost vector 
field. The idea is to compute the renormalized two-point function at one loop approximation 
for the Poincare invariant vacuum state (that obviously verifies (R4) when radiative corrections 
are disregarded) and check whether it satisfies the requirement (R4), while (R1)-(R3) being au- 
tomatically true for a real quantum scalar field. If it is the case, taking ([3| into account with io>2 



given by the one-loop renormalized two-point function, Theorem 1.2 authorizes one to conclude 



that the Hawking radiation viewed as a tunnelling phenomenon survives the introduction of a 
|i</> 3 self-interaction, at one loop at least. 

Before going on, let us briefly discuss why we expect that this simple model makes sense from 
a physical viewpoint in comparison with the analogous situation for a Schwarzschild black hole. 
In the latter case, dealing with the Kruskal manifold, the physically interesting state is the cel- 
ebrated Unruh state, since it is the natural state where Hawking radiation is detected at future 
null infinity. An explicit rigorous construction of that state has recently been established in 
[DMP09] where, in particular, the state has been shown to be of Hadamard type, so it verifies 
the requirement (R4). Moreover, as it was already known from heuristic constructions, that 
state looks like Minkowski vacuum as soon as one approaches the Killing horizon. On the other 
hand the very geometry of Kruskal manifold locally approximates Minkowski one as soon as one 
approaches the Killing horizon. The differences appears far from the Killing horizon. In the 
Schwarzshild manifold the Killing field defining the Killing horizon becomes the Minkowki time 
vector giving rise to the natural notion of the energy far away from the black hole. Instead, in 
Minkowski space, the Killing field defining the Rindler horizon does not approach the Minkowki 
time vector far from the horizon. However, it might not matter since we are interested in what 
happens close to the horizon. Our approximation allows us to exploit the relatively simpler 
version of renormalisation procedure in flat spacetime than the generlly locally covariant version 
in curved spacetime [HoWaOl, HoWa02, BDF09J. A large number of counter terms arising from 
the curvatures will be completely neglected in our elementary model. Certainly, a quantum 
state and the renormalisation procedure are non-local concepts so, there is no guarantee to au- 
tomatically extend a positive result found in Minkowski space for the Poincare invariant state 
to the Kruskal manifold and referring to the Unruh state, although the only local structure 
of the two-point function seems to be relevant. However, if the black hole radiation (viewed 
as a tunnelling probability) did not survive the introduction of a self-interaction in Minkowski 
spacetime, it very hardly would do in curved spacetime. 

2.2. Computation. In the following <j) and <po respectively denote the renormalized and free 
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(massive Klein-Gordon) quantum field. The same convention is exploited for Minkowski vacua, 
and respectively. The Gell-Mann Low formula for time-ordered two-point functions holds: 

where formally: 

S(9)=I + i ! ^^:(n)rf 4 n-i / / : ^ :(tl):0 g :(tl / )(i * fld V + • • • . (7) 

Above we assume that the Wick monomials, denoted by : • : are those defined with respect to 
(their expectation values vanish on \I/o) and the function g £ Cg°(M), attaining constantly 
the value go in a bounded region, has to be switched to an everywhere constant function at the 
end of computations. This is done to remove infrared divergences. Later, we will extract the 
two-point function from the time-ordered two-point function. However this would not be truly 
necessary for, if x does not belong to the causal past of y: 

oj 2 (x, y) := (j){x)<j>{y)^) = (*, T[<t>{x)<f>{y)]^) . (8) 

This is the very situation when x and y are separated by a Killing horizon and x stays in the 
non-static region while y stays in the static one (the right Rindler wedge in our case). 
In the following, we will compute the the explicit expression of the one-loop renormalized two- 
point function in the position domain, rather than in the momentum space, because we intend 
to check whether or not the requirement (R4) still holds taking the radiative corrections into 
account. Here, we are looking for a very precise expression avoiding formal computations based, 
for instance, on divergent series or large momentum approximations. It is of course already 
known that the singularities of two-point function for the theory with radiative corrections are 
for points which are light-like related. However, since such a distribution does not solve the 
Klein Gordon equation, the knowledge of the form of its wave front set is not sufficient to apply 
the result of Radzikowski [Ra96j to conclude that its small distance behavior is of the form (R4). 

Looking at the right hand side of Q one realizes that the first non-trivial contribution to 
the right-hand side of Q comes from the double integral, that is the fish diagram, because the 
previous term yields a vanishing contribution. (In curved spacetime adopting the general locally 
covariant notion of Wick polynomial of [HoWaOl] or even in Minkowski spacetime referring the 
normal ordering to a state different form Minkowski one, also the second term in the right-hand 
side of Q would give a non-trivial contribution.) The relevant 4-point function is therefore the 
distribution on C£°(M 4 ) corresponding to the fish diagram: 

T(x,y,u,u') := (*o,r[^,(x)^o(y):^:(«):^:(«')]*o) • (9) 

As is well-known this is not a well-behaved distribution since it is well-defined only for test 
functions whose supports do not intersect the diagonals of the product M x M x M x M. 
The extension procedure and the classification of the arising ambiguities is nothing but the 
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Epstein- Glaser version of the renormalisation procedure [E]373]> that has been generalized in 
curved spacetime to a generally locally covariant procedure |BF00} IHoWaOl| IHoWa02[ IBDF09] . 
The ambiguities, i.e. the finite-renormalization counter terms, are classified imposing constraints 
concerning, covariance, causal factorisation, scaling behaviour and polynomial dependence on the 
mass and the inverse metric. Dropping terms vanishing in the adiabatic limit, and remembering 
that x 7^ y by hypotheses, the only counter term for the above t(x, y, u, vf) is proportional to 

5T(x,y,u,u) := 5(u,u'){^ ,T[:(pl:(u)(j) (x)(p (y)}^ } . 

If Gp = i(*f>o, r[(/>o(x)0o(y)] l I / o) is the free Feynman propagator we therefore have: 

(¥,T[0(x)0(y)]¥) = -iG F (x,y) 



A g(u) GF(x,u)Gp(u,y)d u 

J M 



2! 



g{u)g{v!) G F {x, u)G F (ext \u, u)G F (z', y)d A ud A u' + ■■■ (10) 



M 2 



where A is a finite renormalisation constant. The first convolution in right hand side is well 
defined, as it follows by direct inspection in flat spacetime or in view of general theorems on 
microlocal analysis in curved spacetime [HoWaOli HoWa02]). G 2 F is well defined as a distribution 
only on C^°(MxM\A2), where A2 := {(x, x)\x G M}, again it follows either by direct inspection 
in flat spacetime or in view of general theorems on microlocal analysis in curved spacetime. Thus, 
in the second convolution in (loj), it has been replaced for an extension G 2 p £Xt ^ acting on the 
whole C^°(M x M). We throughout use the conventions as in [StQ9] about Fourier transform 
(f(x) = (2ir)~ 2 f e*^' 1 /(A;)d 4 A;) and propagators and decompose four vectors as s = (s°,s) with 
s 2 := —(p ) 2 + s • s. Passing to the Fourier transform of distributions, a well-known extension 
of G 2 F is given by: 



G 



2{ext) 



(k) 



2tt 



p 2 + m 2 



ie (p + k) 2 + m 2 



if 



(p 2 + m 2 



d 4 p 



Above and henceforth, a final weak limit indicated by the standard e perscription appears 
explicitely. Making use of the standard Feynman parameters procedure, after a lengthy integral 
computation we obtain: 



G 



2(ext) 



(k) 



2(2ttY 



-1 + A/1 + 



Am 2 



ie 



k 2 



coth 



-1 



1 + 



Am 2 



ie 



k 2 



Inserting the result in (10), exploiting the convolution theorem and performing the adiabatic 



limit, since everting is well defined, so that g is constant, we find: 



-1 



I g 2 GF(x,u)G 2 p ext \u,u')GF{u ' : y)d A ud i U 
2! J M 2 



4„./ 



4(2vr) 



ikfi(x—y)i* 



{k 2 + m 2 



+ A 1 + 



Am 2 



ie 



k 2 



coth 



-1 



1 + 



Am 2 



ie 



k 2 



d A k 



(11) 
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k° 



Figure 1: The picture emphasizes the two poles at z = ±P = ±Vk 2 + m 2 — ie and the two 
branch cuts (starting at z = ±B = ±\/k 2 + 4m 2 — ie) relevant for the complex extension z of 
k° G R. It represents also the two contours on which the 2 integration is taken for positive 
(lower contour) and negative (upper contour) t x — t y . 



The integration in k° can be computed before that in k extending k° to a complex variable z. 
It is done by completing the integration along the real line into a closed contour with an arch 
at infinity in either the lower or the upper half-plane, depending on the sign of (t x — t y ), taking 
advantage of the residue technology as well known. The analytic continuation (in the variable k°) 



of the integrand of ( 11 ) gives rise to a couple of poles at ±\/k 2 + m 2 — ie respectively. However, 



although no problems arise with the decay rapidity of the considered functions on the portion 



of the contour at infinity, a closer scrutiny of the function in square brackets in (11) added to 
— 1 reveals the appearance of a pair of logarithmic branch cuts. One, relevant for t x — t y > 0, 
completely stays in the lower half-plane starting from \/k 2 + 4m 2 — ie and joining oo; the other, 
relevant for t x — t y < 0, completely stays in the upper half plane starting from — Vk 2 + 4m 2 — ie 
and joining oo. So the contributions of these branch cuts have to be taken into account. Making 
explicit the contribution of the poles, re-adjusting the e-prescription into an equivalent form, 
and indicating by C(x,y) the contribution due of the branch cuts, the result is: 

2 f 

- — ' G F{x,u)G 2 p ext \u,u)G f{u \y)d 4 ud i u 



2! 

9 2 



M 2 



4(2vr) 



K mJ a(x, y) + 2i\t x - t y \e + e 2 



2 /i \ K\ (mJa(x,y) + 2i\t x - t y \e + e 2 

-ttL i-^F V ; ^ + C(x,y) (12) 

4(2vr) 2 V2 3^3/ my/a(x, y) + 2i\t x - t y \e + e 2 V ' 

where K v are the well-known modified Bessel functions of the second kind. We have used 
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formulas 3.914 (10) and (9) of |GR07] . The term proportional to the undetermined constant A 
in ( 10 ) can be analogously computed and the only result is to change the factor in front of Kq by 



an unknown constant. The term in (|12[) immediately before C(x,y), if replacing the coefficient 



in front of K\ by i ^ is nothing but Gjr(x,y). Putting all together, (10) yields the result: 



T[0(x)<Ky)]tf > = AK o mJa(x, y) + 2i\t x - Ue + e 



+ 



m 



(2tt) 2 4(2vr) 



— ) 
3^3/ 



Ki (my/a(x, y) + 2i\t x 



Ue + e* 



m^fo{x,y) + 2i\t a 



t v k + e 2 



+ C(x,y) + 



(13) 



It remains to evaluate C(x, y). Using the very definition of coth 1 and the well-known fact that 
ln(— \x\ +iS) — ln(— |x| —iS) — > 2m for 5 — > + in evaluating the integral along curve surrounding 



a branch cut of the function added to —1 in square brackets in (11), one finds: 



C(x,y) 



4(2vr) 



d|k||k| 



dz- 



(z 2 - k 2 



sin(|x 



m 2 + ie) 2 



1 + 



Am 2 



ie 

~7i 



where 7 is the lower branch cut. However, since the integrand of the z integration is holomorphic 
in the lower half plane (barring a branch cut from z = |k| to z = \/k 2 + 4m 2 — ie) and it decreases 
rapidly, 7 can be deformed without affecting the value of the integral, provided the path keeps 
joining \/k 2 + 4m 2 — ie and 00 (and avoides the cut). For convenience we therefore assume 7 
to be of the form z(s) = \/k 2 + 4m 2 — ie + s 2 with s G [0, +00). The integration in ds can be 
evaluated after computing that in d|k| obtaining: 



C(x, y)=X (a(x, y) + 2i\t x - t y \e + e 2 ) , 
where M 2 = s 2 + 4m 2 and: 



%{u) 



9" 
4(2vr) 4 



2m 



4m 2 



(M 2 



rn 



2\2 



MKi (MVm) 



(14) 



(15) 



We stress a feature shared by all the three functions of x and y in the right-hand side of (13) 



taking (14) and (15) into account. When x and y are spacelike related, a(x,y) > so that the 
part 2i\t x — t y \e + ? of a(x, y) + 2i\t x — t y \e + e 2 does not affect the final result when taking the 



weak limit e — > . When x stays in the (causal) future of y, \t x 
x in the internal region and y in the external region - so that 



t 



y\ °x 
holds true 



t y . Therefore, for 
- we can replace 



a(x,y) + 2i\t x — t y \e + e for the regularized distance appearing in the Hadamard prescription 



(7, 



(x, y) = a(x, y) + 2i(t x - t y )e + e 



1 We can thus write, with x and y as stated above: 
OM(aMy)tf) = 



10 



AK Q L/m 2 cr £ (x,y) 



m 



(2vr) 2 4(2^) 2 



1 7T \-\K 1 ^m 2 a e (x,y) 

2 3\/3/J y / m 2 a e (x,y) 



+ X{a e (x,y)) + ...(16) 



Taking the complex conjugate of both sides of ( |16| ), using the fact that <f> is Hermitian, and fi- 
nally interchanging the name of x and y, from the elementary properties of the K u functions, one 



easily see that ( 16 ) holds also for x in the causal past of y. So (16) holds for all values of x and y 



2.3. Hawking radiation in the tunnelling approach survives the interaction. We 



intend to analyse the short-distance behaviour of the right-hand side of ( 16 ) to check if it fits the 



requirements (R4), especially taking Remark |1.1| into account. Let us start by considering the 
last term in the right-hand side of ( |16[ ) which deserves more attention. First of all we recall to the 
reader that the function K%(Q — 1/C is bounded in the closed half-plane ReC, > 0. Boundedness 
away from the origin easily follows from 8.451(4) and 8.451(6) in [GR07J, while boundedness 
around the origin is consequence of the decomposition: 

Kx{Q = i + /i(C)ln(C/2) + ^i(C) CeC (17) 

where the modified Bessel function I\ and tpi are holomorphic in the whole complex plane with 
li(0) = 0. The definition of X in ([15} yields: 



X(a e 



D h e 

Of ^/ol 



where 



h t (x,y) 



D 



9 



4(2vr) 4 

„2 



9 



4(2?r) 4 



+ " VP" 



4m 2 



m 



2\2 



M 



1 



4m 2 



m 



2\2 



4(2vr) 2 m 2 



(18) 

(19) 
(20) 



Since Re (^/o T l) > 0, the function in squared brackets in the integral in the right-hand side of 
(19) is bounded. Consequently h e satisfies the hypotheses stated in (1) in Remark |1.1| The first 
term in the right-hand side of (16) can be treated exploiting the known expansion: 



^o(C) = -/o(C)ln(C/2) + Vo(C) CeC 

where the modified Bessel function Iq and ipo are holomorphic in the whole complex plane. The 



It is worth noticing that, due value of D in (20), the leading divergence 
in (18) exactly cancels an analogous divergence proportional to g 2 arising form the 



second term in the right-hand side of (16) can analogously be treated taking advantage of (17). 
Everywhere £ 
of %{a, 



second term in the right-hand side of (16). Collecting all the contributions together, we can 
easily conclude that 



uj 2 (x,y) 



1 

(270Ve 



+ 



h e (x,y) 
V a e{x,y) 



+ fe(x,y)lna e (x,y) +r e (x,y) , 



(21) 
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where the functions h € ,f £ ,r e verify the conditions stated in (1) in Remark Regarding f e 
and r e , it follows from the fact that the Taylor expansions of I\ and ipi centered on the origin 
are made of odd powers of £ only, while those of Iq and ipQ are made of even powers of £ only. 



We stress that the leading term in the right-hand side of (21) is the same as in the free theory 
because, as we have found, the radiative correction give no contribution to the dominant diver- 
gence of the two-point function. 

Since L02 satisfies the requirement (R4), we can apply Theorem |1.1| obtaining that: 

Even taking the radiative corrections of the interaction Lj = ||</> 3 into account, i.e. referring to 
the renormalized vacuum state \T/ and renormalized field operators eft, at one loop: 

linj l|<*M(/a)<K/a)*>| 2 ~ CEl e-P* E ° , 

A— >0 

for packets sharply concentrated around a large value Eq of the energy. The factor C includes 
contributions due to the self-interaction. 



1.1 



We remark that also the second case in Theorem 
when both / and /' have supports in the external region. 



holds for K¥,0(/aM/a)¥}| 2 , that is (1.1) 
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